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An effect of rotation on a developed turbulent stratified convection is studied. Dependencies of 
the hydrodynamic helicity, the alpha-tensor and the effective drift velocity of the mean magnetic 
field on the rate of rotation and an anisotropy of turbulent convection are found. It is shown that 
in an anisotropic turbulent convection the alpha-effect can change its sign depending on the rate of 
rotation. The evolution of the alpha-effect is much more complicated than that of the hydrodynamic 
helicity in an anisotropic turbulent convection of a rotating fluid. Different properties of the effective 
drift velocity of the mean magnetic field in a rotating turbulent convection are found: (i) a poloidal 
effective drift velocity can be diamagnetic or paramagnetic depending on the rate of rotation; (ii) 
there is a difference in the effective drift velocities for the toroidal and poloidal magnetic fields; (iii) 
a toroidal effective drift velocity can play a role of an additional differential rotation. The above 
effects and an effect of a nonzero divergence of the effective drift velocity of the toroidal magnetic 
field on a magnetic dynamo in a developed turbulent stratified convection of a rotating fluid are 
studied. Astrophysical applications of the obtained results are discussed. 

PACS numbers: 47.65.-|-a; 47.27.-i 



I. INTRODUCTION 

Turbulent transport of particles and magnetic fields 
was intensively studied for the Navier-Stokes turbulence 
(see, e.g., ||, ^, ||). However, there are a number of 
applications with other kinds of turbulence, e.g., turbu- 
lent convection. For instance, in the Sun and stars there 
is a developed turbulent convection that is strongly in- 
fluenced by a fluid rotation. 

The mean-field theory of magnetic field was in gen- 
eral developed for the Navier-Stokes turbulence with- 
out taking into account turbulent convection (see, e.g., 
§ i 0> i S lO' 111)- III particular, the dependen- 
cies of the the a-effect, the effective drift velocity and 
the turbulent magnetic diffusion on the rate of rotation 
were found only for the Navier-Stokes turbulence (see, 
e.g., |lj, |l^) in spite of that in many astro- 

physical applications there are turbulent convection re- 
gions. A turbulent convection in different situations has 
been studied mainly by numerical simulations (see, e.g., 
0,§|,^,|^,|l|,||,||). 

In this paper we study an influence of rotation on a 
developed turbulent stratified convection. This allows us 
to find the dependencies of the hydrodynamic helicity, 
the alpha-tensor and the effective drift velocity of the 
mean magnetic field on the rate of rotation. 

This study has a number of applications in astro- 
physics. In particular, the evolution of the mean mag- 



netic field in the kinematic approximation (without tak- 
ing into account a two-way coupling of mean magnetic 
field and turbulent fluid flow) can be described in terms 
of propagating waves with the growing amplitude, i.e., 
the magnitude of the mean magnetic field B is given by 
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where Bq is a seed magnetic field, 7^ is the growth rate of 
the mean magnetic field, tOg and k are the frequency and 
the wave vector of a dynamo wave. In the Sun, e.g., ac- 
cording to the magnetic field observations these dynamo 
waves with the ~ 22 years period propagate to the equa- 
tor (see, e.g., ||, ^, 0, |l^). The magnetic field is gener- 
ated in the turbulent convective zone inside the Sun. The 
growth of the mean magnetic field is a combined effect 
of a nonuniform fluid rotation (the differential rotation, 
V((5ri)) and helical turbulent motions (the a-effect). The 
direction of propagation of the dynamo waves is deter- 
mined by a sign of the parameter a[d{SQ)/dr], where 
r, 9, if are the spherical coordinates, and is the angu- 
lar velocity. When the parameter a [d{6fl) /dr] is nega- 
tive the dynamo waves propagates to the equator. The 
helioseismology shows that in the solar convective zone 
d{Sn)/dr > and the existing theories yield a > 0. This 
results in that the dynamo waves should propagate to the 
pole in contradiction to the solar magnetic field observa- 
tions (see, e.g., [| | 0, |, |ig). 

In this study we found that in a developed turbulent 
convection the a-effcct can change its sign depending on 
the rate of rotation and an anisotropy of turbulence. In 
the lower part of the solar convective zone the fluid ro- 
tation is very fast in comparison with the turnover time 
of turbulent eddies. In this region a > 0. In the up- 
per part of the solar convective zone the fluid rotation is 
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very slow and a < 0. This explains the observed proper- 
ties of the solar dynamo waves. The growth of the mean 
magnetic field is saturated by nonlinear effects (see, e.g., 
p^ , I25I , ^ ^). The 22-years solar magnetic activ- 
ity is also poorly understood. A characteristic time of 
the turbulent magnetic diffusion in the solar convective 
zone is of the order 2-3 years and it cannot explain the 
characteristic time of solar magnetic activity. We found 
that the fast rotation causes an additional effective drift 
velocity of a mean magnetic field that can increase the 
period of the dynamo waves provides the 22-years solar 
magnetic activity. 



II. THE GOVERNING EQUATIONS AND THE 
METHOD OF THE DERIVATIONS 

Our goal is to study an effect of rotation on a developed 
turbulent stratified convection. This allows us to derive 
dependencies of the hydrodynamic helicity, the alpha- 
effect and the effective drift velocity of the mean magnetic 
field on the angular velocity. To this end we consider a 
fully developed turbulent convection in a stratified ro- 
tating fluid with large Rayleigh and Reynolds numbers. 
The governing equations are given by 
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where u is the fluid velocity with V-u = A- u, D/Dt ~ 
d/dt -t- u • V, 11 is the angular velocity, g is the gravity 
field that includes an effect of the centrifugal force, p^i^, 
is the viscous force, F^{S) is the thermal flux that is 
associated with the molecular thermal conductivity, A — 
-p-iVp„, Nb = (7Po)"^VPo-Po"^Vpo. The variables 
with the subscript "0" corresponds to the hydrostatic 
equilibrium (i.e., the hydrostatic basic reference state): 



VPo =PoS 



(4) 



and To is the equilibrium fluid temperature, S = P/^Pq — 
p/Po are the deviations of the entropy from the hydro- 
static equilibrium, P and p are the deviations of the fluid 
pressure and density from the hydrostatic equilibrium. 
The Brunt- Vaisala frequency, Jib, is determined by the 
equation Ql — — g • N;,. To derive Eq. (^) we use an 
identity: -VP + gp = -po[V{P/po) + gS - PN^/po], 
where we assumed that |PN{,/po| < \gS\, iPN^/pol < 
|V(P/po)|- This assumption corresponds to nearly isen- 
tropic basic reference state when Nf, is very small. For 
the derivation of this identity we also used Eq. (^. We 
also consider a low-Mach-number fluid flow with a very 
small frequency fib, i.e., \Qb\ <C \^g^ and |f2bTp <C 1, 
where r is the correlation time of the turbulent velocity 
field. Equations (|^) and (^) are written in the Boussinesq 
approximation for V • u 7^ 0. This is more usually called 
" the anelastic approximation" . 



Now we consider a purely hydrostatic isentropic basic 
reference state, i.e., O;, = (N^ = 0). Thus the turbu- 
lent convection is regarded as a small deviation from a 
well-mixed adiabatic state (for more discussion, see pof). 
We will use a mean field approach whereby the veloc- 
ity, pressure and entropy are separated into the mean 
and fluctuating parts. Using Eqs. ^ and (^) we de- 
rive equations for the turbulent flelds: Vz = \/ po{z) Uz, 
w = a/ po{z) (Vxu)2 and s = ^ pq{z) {S — S), where 
S = (S) is the mean entropy, the angular brackets de- 
note ensemble averaging, and for simplicity we consider 
turbulent flow with zero mean velocity. Here po{z) is di- 
mensionless density measured in the units of po{z ~ 0). 
The equations for the turbulent fields are given by 
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where VL^ 



-g -VS, Aj^ = A - d^/dz^, Vn, Wn and 
Sn are the nonlinear terms which include the molecu- 
lar dissipative terms [see Eqs. ( |A10| )-( |A12|) in Appendix 
A], the field g is directed opposite to the axis z and 
n = (Slj,, 0, f^z). We assumed here that A" 1 1 9A/92: 1 < A. 
Equation (p) follows from Eq. (^) after the calculation 
[Vx(Vxu)]2. 

By means of Eqs. (^-(0) we derive dependencies of the 
hydrodynamic helicity, the alpha-effect and the effective 
drift velocity on the angular velocity. The procedure of 
the derivation is outlined in the following (for details, see 
Appendixes A, B and C). 

(a). Using Eqs. (|^)-(0) we derive equations for the 
following second moments: 

/jj(k) = L{v.i,Vj) , x(k) = L{w,Vz) , 



P(k) = Lis,w) , G(k) 
$,(k) = L{s,v,) , e(k) 



L{w,w) , 
L{s,s) , 



where L{a,b) = (a(k)&(— k)) and v — ^ Pq{z) u. The 
equations for these correlation functions are given by Eqs. 
(|A4|)-( |A9| ) in Appendix A. In this derivation we assumed 
that A^ P. 

(b). The equations for the second moments contain 
third moments and a problem of closing the equations for 
the higher moments arises. Various approximate meth- 
ods have been proposed for the solution of problems of 
this type (see, e.g., [|l|, |3l|). The simplest proce- 
dure is the T approximation, which is widely used in the 
theory of kinetic equations. For magnetohydrodynamic 
turbulence this approximation was used in ]3^ (see also 
p7| , ^ HI). One of the simplest procedures which al- 
lows us to express the third moments /^^ , Xn j ■ ■ • 1 in 
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Eqs. 



)-(|A9D in terms of the second moments, reads 
/(k)-/("'(k) 



/„(k)-/f (k) 



T{k) 



(8) 



and similarly for other third moments, where /(k) = 
BiCj fijCk), e is the unit vector directed along the axis z, 
the superscript (0) corresponds to the background turbu- 
lent convection (it is a turbulent convection without ro- 
tation, r2 = 0), and T{k) is the characteristic relaxation 
time of the statistical moments. Note that we applied 
the T- approximation only to study the deviations from 
the background turbulent convection which is caused by 
the rotation. The background turbulent convection is 
assumed to be known. 

The T-approximation is in general similar to Eddy 
Damped Quasi Normal Markowian (EDQNM) approx- 
imation. However, there is a principle difference between 
these two approaches (see @ The EDQNM clo- 

sures do not relax to equilibrium, and this procedure does 
not describe properly the motions in the equilibrium state 
in contrast to the r-approximation. Within the EDQNM 
theory, there is no dynamically determined relaxation 
time, and no slightly perturbed steady state can be ap- 
proached In the r-approximation, the relaxation 
time for small departures from equilibrium is determined 
by the random motions in the equilibrium state, but not 
by the departure from equilibrium [ pO[ . As follows from 
the analysis performed in [30| the r-approximation de- 
scribes the relaxation to the equilibrium state (the back- 
ground turbulent convection) more accurately than the 
EDQNM approach. 

Note that we analyzed the applicability of the r- 
approximation for description of the mean-field dynam- 
ics of the mean magnetic field and mean scalar fields 
by comparison of the derived mean-field equations using 
other methods such as the path-integral approach and the 
renormalization group approach (see ||,|36[|7[|38[ |9). 
This comparison showed that the r-approximation yields 
the results similar to that obtained by means of the other 
methods. 

(c) . We assume that the characteristic times of vari- 
ation of the second moments /(k), x(k),..., 0(k) are 
substantially larger than the correlation time r(A;) for all 
turbulence scales. This allows us to determine a station- 
ary solutio n for t he sec ond moments /(k), x(k), . . . , 9(k) 
[see Eqs. ([A2^)-( |A30D in Appendix A]. 

(d) . For the integration in k-space of the second mo- 
ments /(k), x(k), . . . , 0(k) we have to specify a model 
for the background turbulent convection (without rota- 
tion) . Here we use the following model of the background 
turbulent convection which will be discussed in more de- 
tails in Appendix D: 



$r^(k) 



+ii^(o)(k)(exk),] , 
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$i°)(k) = <p:w{k) 



F(")(k) 
G(°)(k) 

e(")(k) 



-6*(**-(exk))/W(k)//, , 
(1 + £)/(") (k)fc2 ^ 

2e^W{k) , 



(11) 
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where /y(k) = (v,{k)vj{-k)), W{k) = W{k)/8Trk^ 
/(0)(k) = Mk^/k^Wik), and /^)(k)e., = /(o)(k). 



is the degree of anisotropy of the turbulent velocity field 
u = + u^e. Here Pij{k) — Sij — kij, kij = kikj/k'^, 



k = k^ + fc,e, fc, = k • e, Pl^\k^) = 6, 



(k_L),(ki)j/fci, e. 



k 



BiBj, a is the degree of 



anisotropy of the turbulent flux of entropy (see below). 
We assume that r(fc) = 2r„f(fc), W{k) = —df{k)/dk, 
f{k) = (fc/fco)"'^"*, 1 < (7 < 3 is the exponent of the 
kinetic energy spectrum (e.g., q — 5/3 for Kolmogorov 
spectrum), ko = l/^o, and is the maximum scale 
of turbulent motions, r^ = Iq/uq and uq is the char- 
acteristic turbulent velocity in the scale Iq. Motion in 
the background turbulent convection is assumed to be 
non- hehcal. In Eqs. @ and 1^) we neglected small 
terms ^ 0(A/*) and ^ 0(A<I>*), respectively. Now 



we calculate f^j^ = J f^^' (k) dk using Eq. 



(0), 



^(0) ^ 

J ij 

{f*/3)[dij + (3e/4)((5y - cy)]- Note that -4/3 < £ < oo. 
The lower limit of s follows from the condition fx^J > 
(or > 0). Similarly, using Eqs. (0)-(|2|) we obtain 
$(0) = J 4>(")(k) dk = The parameter a can be pre- 
sented in the form 



l + |(g + l)/(g-l) 

1 + e/s 

e = iijhy-' - 1 , 



(15) 
(16) 



where l± and Iz are the horizontal and vertical scales 
in which the two-point correlation function ^f "^ (r) = 
(s(x)u(x -I- r)) tends to zero. The parameter ^ deter- 
mines the degree of thermal anisotropy. In particular, 
when l± — Iz the parameter ^ = and a = 1. For ^ Iz 
the parameter ^ = —1 and a — — 3/(g— 1). The maximum 
value ^niax of the parameter ^ is given by ^max = 9—1 for 
CT = 3. Thus, for ct < 1 the thermal structures have the 
form of column or thermal jets (Z_l < Iz), and for a > 1 
there exist the "pancake" thermal structures {l±^ > Iz) va. 
the background turbulent convection. 

T he r elationship between $* and /, follows from 
Eq. ( |Al| ) for the kinetic turbulent energy p^^{\}?), and it 
is given by /, = 2XgTg^l/e, where A = 2eSr/ (e + 2) and 
(5, = (3 — q)/'^{q — !)■ Note that for Kolmogorov spec- 
trum (7 = 5/3 and (5* = 1. In Section HI we will present 
results for J* = 1. For the integration in k-space we used 
identities given in Appendixes B and C. 
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Thus, the "input parameters" in the theory include the 
parameters that describe the model of background turbu- 
lent convection, i.e., the degree of anisotropy of the tur- 
bulent velocity field e, the degree of anisotropy of the tur- 
bulent flux of entropy a, the maximum scale of turbulent 



and for w ;:!> 1 it is given by 



motions Iq, the turbulent velocity uq = -^Z (u^) — 
(the r.m.s. velocity) in the maximum scale of turbulent 
motions, the exponent of the kinetic energy spectrum q. 
The "input parameters" also include the density strati- 
fication length A and the angular velocity $1. Note that 

= ul/{l + e/2) and $* = Mg/(2^*gTj. The described 
above procedure allows us to determine the dependencies 
of the hydrodynamic helicity, the alpha-effect and the ef- 
fective drift velocity of the mean magnetic field on the 
rate of rotation. 

The considered model of a background turbulent con- 
vection written in k-space is enough general and it does 
not contradict to the known Nusselt number dependen- 
cies on Rayleigh number. On the other hand, the ob- 
servations of the turbulent convection on the surface of 
the Sun cannot give the Nusselt number dependence on 
Rayleigh number, i.e., it is possible to obtain only one 
point in this curve. The parameters e, uq, Iq, $7, etc can 
be calculated from the solar observations. In addition, 
the direct numerical simulations of turbulent convection 
(see HI) are in an agreement with our model of 

turbulent convection. 



III. EFFECT OF ROTATION 

In this Section we present the results of the calcula- 
tions (described above) for the hydrodynamic helicity, 
the alpha-effect and the effective drift velocity of the 
mean magnetic field as the functions of the rate of ro- 
tation and an anisotropy of turbulence. 



A. The hydrodynamic helicity 



Using Eqs. ( A34| ) and (A3J) in Appendix A we find 
the dependence of the hydrodynamic helicity x*^"-* — 
(u • (Vxu)) on the angular velocity: 



1 

12 



[^!i{lo) + ^2{i^)siv? (f>i 



-f'^'3(cj)sin <j)i] sin <j)i 

(for details, see Appendix A), where u 
MoT,,, Uq — 2(7Tg$*(5, , sin0i — u} ■ e, is the latitude, 
u; = n/n, e is the unit vector directed along the 2:-axis 



la 



Lp = A ^, the functions are given by Eqs. (CI) 



in Appendix C. Hereafter we assume that (5, = 1. For a 
slow rotation {oj ^ 1) the hydrodynamic helicity x^'") is 
given by 
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(19) 



Note that the meaning lu = 4nTj, 1 is w large, but 
only up to some upper limit, i.e., an intermediate range 
of values. This implies that the rotation cannot be very 
fast to affect the correlation time T(fc) of turbulent veloc- 
ity field in its inertial range. Also we assumed that the 
parameters e and a are independent of lo. 



B. The a-effect 

Now we find the dependence of the a-effect on the an- 
gular velocity. To this end we use the induction equation 
for the magnetic field 



an 

'dt 



V X (u X H - 77V X H) 



(20) 



where i] is the magnetic diffusion due to the electrical 
conductivity of fluid. The magnetic field, H, is divided 
into the mean and fluctuating parts: H = B -I- b, where 
the mean magnetic field B = (H) and b is the fluctuating 
field. An equation for h = -/pob follows from Eq. ([2^ 
and it is given by 



dh 

'dt 



(B • V)v - (v • V)B - (v • A)B 



-(B- A)v + H 



N 



(21) 



where H^r are the nonlinear terms which also include 
the magnetic diffusion term [see Eq. (A47) in Ap- 
pendix A]. In order to derive equation for the a-tensor 
we introduce the electromotive force £i = (u x b)^ = 
Po^Sirnn J Xmn(k)dk, whcrc X^if 0^) = ('^i (k)^j ("k)) = 
L{vi,hj). A general form of the electromotive force 
is given by £i = atjBj + (V'=^xB)i - r]ij{VxB)j - 
kidB)ij - [<5x(VxB)]i = aijBj + bijkBij (see, e.g., 
and Appendix A), where the tensors a^- and rjij de- 
scribe the a-effect and turbulent magnetic diffusion, re- 
spectively, is the effective diamagnetic (or paramag- 
netic) velocity, Kijk and S describe a nontrivial behavior 
of the mean magnetic field in an anisotropic turbulence, 
B,,, = V,B, and (9B)y = (l/2)(B,j + 5^-,). The a- 
(17) tensor, a^, is determined by a symmetric part of the 



tensor a^, i.e., by a\j = {l/2){aij +aji). The tensor Uij 
is calculated in Appendix A. The a-tensor is given by 

1 

a,j = -i^-j- j {sin (1)1 [{-^iiuj) + '<if5{uj) sin^ (l)i)S,j 

+ («'6(t^) + 4'7(w) sin^ (t)i)LJij + *8(t^) e,j] 

+ [1'9(w) + ^w{lu) sin^ </>/] (e.Ldj + ejW,)} , (22) 
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(for details, see Appendix A), where uj. 



<ij — OJiUJj, 



BiBj, the functions ^'^(ti;) are given by Eqs. ( pl| ) i 
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FIG. 1: Characteristic ranges of parameters with differ- 
ent behavior of a-effect (FIG. la) and the parameter ~ 
-(l/3)r(,x'''' (FIG.lb). The range I for the a-effect (FIG.la) 
also exists for —2 < A < and —9/2 < a < 3. 

Appendix C. Here we present asymptotic formulas for 
the isotropic part = aSij) of the a-tensor. For a 

slow rotation (w <C 1) the parameter a is given by 

and for w ^ 1 it is given by 

+ {(T -I) SIT? 0; j sin (j)i . (24) 

It is seen from Eqs. ( p^ ) and ( p3| ) that for a slow rotation 
and isotropic background turbulent convection [a — 1 
and e = 0), the parameter a « — (5/27)T(,x^"\ where 
X^"-' — (u • (Vxu)). However, when a rotation is not 
slow, the latter relationship does not valid. 

The Qf-effect depends on the degrees of the velocity 
anisotropy e and the thermal anisotropy cr. Asymptotic 



formulas for a slow rotation (w <C 1) and for a; ^ 1 
show that there are several characteristic ranges of pa- 
rameters with different behavior of a-effect. In FIG. la 
these ranges are separated by lines a = 3(3 — 2A), 
cr = 3(1 - A/2) and o- = 6 - 5A/2, where -2 < A < 2, 
-9/2 < cr < 3, and < 0; < 7r/2. Here A = 2£/(e + 2). In 
the ranges I and II the a-effect does not change its sign 
for all r^TQ and (pi. In particular, in the range I: a > 
and in the range II: a < 0. In the range V the a-effect 
changes its sign at a certain value of fir^ for all pi . In the 
ranges HI and IV the a-effect changes its sign at a cer- 
tain value of SItq and a certain range of the latitudes pi . 
In the range HI the degree of thermal anisotropy cr > 1 
(which corresponds to the "pancake" small-scale thermal 
structure of the background turbulent convection), and 
in the range IV the degree of thermal anisotropy cr < 1 
(i.e., a column- like thermal structure). The a-effect can 
be negative for a slow rotation only in the range II. Note 
that the negative a-effect corresponds to the propagation 
of the solar dynamo waves to the equator. 

Our analysis shows that when the rotation is not slow, 
the a-effect is determined not only by the contributions 
from the hydrodynamic helicity and its behavior is much 
more complicated in a rotating fluid. In order to demon- 
strate this we plotted in FIGS. 2-4 the dependencies of 
the a-effect (solid line) and = —{l/3)Tgx''^^ (dashed 
line) on the parameter Ut^ for different latitudes (FIG. 
2 is for the latitude = 15°, FIG. 3 is for pi = 35° 
and FIG. 4 is for pi = 90°). Here the parameters a and 
a^ are measured in the units of laUo/4:Lp. Figures 2-4 
demonstrate that the functions a{flTg) and a^{QTg) are 
totally different. For example, in the case e = 13 and 
(7 = the a-effect and have opposite signs for all Qt^ 
(see FIG. 4c). 

Figure lb shows the ranges of parameters {a and A) 
with different behavior of a^. In FIG. lb these ranges 
are separated by lines a = (9/41)(25A/12 — 1) and A = 0, 
where —2 < A < 2 and —9/2 < cr < 3. The numeration 
of the ranges in FIG. lb for a^ is the same as for the 
parameter a in FIG. la. Comparison of FIGS, la and lb 
shows that the ranges HI and IV (whereby the a-effect 
changes its sign at a certain value of rir^ and a certain 
range of the latitudes pi) do not exist for a^. On the 
other hand, there is a new range (the range VI) in FIG 
lb whereby the sign of a^ changes from negative value for 
a slow rotation to positive value for a; ^ 1. The locations 
of the ranges II and V for a^ are different from that of 
the a-effect. Therefore, the behavior of the parameter 
and the a-effect are different in a rotating fluid. 

The dependencies of the a effect on the latitude pi for 
different values of the degrees of anisotropy e and a, and 
different values of the parameter rir^ are shown in FIG. 5. 
It is seen in FIG. 5b that the a effect changes its sign at 
pi « 20° — 40° for JItq — 5 (this value of corresponds 
to the lower part of the solar convective zone). 

In view of the applications to the astrophysics, the case 
with negative a-effect for pi > Q is most important be- 
cause this provides a propagation of the solar dynamo 
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FIG. 2: The a (solid line) and = ~(1/3)t-oX*"' (dashed 
line) as the functions of the parameter f2r„ for 0; = 15° and 
different values of the degrees of anisotropy: (a), e = and 
cr = 1; (b). e = 1.2 and cr = 2; (c). e = 13 and cr = 0. In FIG. 
2c the a effect is multiplied by 5. 



waves to the equator according to the solar observations 
(see, e.g., |, |, 0, |, [lO| ) . 



C. The effective drift velocity of the mean 
magnetic field 



Now we determine the effective drift velocity 
— {l/2)ekija[f^^ — V^^^ +V^'^^ of the mean magnetic £ 
using Eq. (A69), where a[f^^ = (l/2)(ajj — aji), and 

v(i) 
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FIG. 3: The a (solid line) and = -(l/3)roX*''' (dashed 
line) as the functions of the parameter ilr^ for (pi — 35° and 
different values of the degrees of anisotropy: (a), e = and 
a = 1; (b). e — 1.2 and a — 2; (c). e = 13 and (7 = 0. 



(for details, see Appendix A), where r, 6, ip are the spher- 
ical coordinates, 0/ = Tr/2 — 6, tDxe = cos^iie^p, uj = 
e sin (f)i — eg cos0/, the functions E^.{lu) are given in Ap- 
pendix C. For a slow rotation (w ^ 1) the effective drift 
velocities are given by 



e 1- 



4 / IqUq 
15\ Lp 



4(6 + 2) 

a 
6 



3£ - 1 
'7(7+2) 



V(2) 



5 V L, 



4 fi^n 



sin(2(/); 

cr 5(1 -e) 
" 6 ^ 9{e + 2) 



(27) 

(u,xe) , (28) 



and for w 3> 1 they are given by 



V(i) « ^^(^J^]le(l 



8uj\ Lp J 



[e(l + sin^ 
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FIG. 4: The a (solid line) and = -(l/3)roX^"' (dashed 
line) as the functions of the parameter JItq for (fn = 90° and 
different values of the degrees of anisotropy: (a). £ = and 
a- = 1; (b). e = 1.2 and cr = 2; (c). e = 13 and a = 0. In FIG. 
4b the a effect is multipUed by 2 and in FIG. 4c by 5/2. 
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FIG. 5: The dependence of the a effect on the latitude 0; 
for different values of the degrees of anisotropy: e = 1.2 and 
(T = 2 (dashed); e = 13 and o" = (solid); e = 13 and a — 2.2 
(dashed-dotted); e = 13 and a = 0.415 (dotted); and for 



different values of the parameter Qr^ 



0.1; (b). 



ilTg = 5. The dashed-dotted line in FIG. 5b shows a/5. The 
latitude is measured in degrees. 



-ee sm{24)i)] 



(29) 



v(^) - i(T7)("-')'^"''^'('^^^)- ^^^^ 

For uj — this effective drift velocity, V*^^), corresponds 
to the well-known turbulent diamagnetic velocity (see, 
e.g., P, ^, 1^, ^, ^). Indeed, since we suggested that 

V(p(lP)) « 0, thus V(u2)/(u2\ - „„A T7„ 

for = reads 



L-^e and Eq. (|2 



4 /'iln 



15 V L 



1 



(7 

6 



4(6 + 2) 



,(31) 



where = Iq/uq. Figure 6 shows the effective drift ve- 
locities: vj^^ and Vy'' as the functions of the parameter 
V,Tg for different values of the degrees of anisotropy. 

The effective drift velocity V*^^' causes an additional 
differential rotation. Indeed, let us introduce the an- 



gular velocity difference, 5fl, which is determined from 
the identity: V*^^^ — Sil,r (d>xe). Comparison of this 
definition with Eqs. (28) and ( pO[ ) yields equations for 
Sil{r) oc r~^. Calculating the r-derivatives of Sn{r) we 
obtain equations which determine the differential rota- 
tion for a slow rotation (lu <ti I) 



d{sn) _ 

dr 

and for w 3> 1: 
dr 



I - - 



5(1 - e) 
9(e + 2) 



1 

2^ 



(1 



(32) 



(33) 



(c) 

aij Bj which for 



The electromotive force has a term 
an axisymmetric case contributes only to an additional 
effective drift velocity, V^^^ , of the mean magnetic field, 
i.e., 



Ac) 



= [V(3)x(Bp-Bt)]. 



(34) 



8 



-18(cr- 1)C0S2 (j>l]} 



(36) 





FIG. 6: The effective drift velocities: (a), vj^' for <j>i = 45°; 
(b). Vr^' for (f)i = 90° as the functions of the parameter Qt^ 
for different values of the degrees of anisotropy: e = 13 and 
a = —2.2 (solid); e = 1.2 and a — 2 (dashed); £ = 13 and 
CT = (dotted); £ = and a = 1 (dashed). 



(for details, see Appendix A), where B = Bt + Bp is 
the mean magnetic field with the toroidal, By, and the 

(c) 

poloidal, Bp, components, the tensor a], is determined 



by Eq. (A59), and the additional effective drift velocity 
is given by 



V(3) 



1 

24 



{e[Ej{uj) + Esiuj) sin^ 
+ l-e0[Eio{Lj) 



+Eq{u}) sin'' 



-Eg{uj) sin^ (?!);] sin(2(?!)i)} 



(35) 



(for details, see Appendix A). Note that for a slow rota- 
tion (w <C 1) the additional effective drift velocity is very 
small, i.e., T^*^^' ^ 0{uj'^)^ and for w 3> 1 it is given by 

V(3) ^ ^fi2!^^cos,/)z{ecos0i[A + lO-13a 
— 18((T — 1) sin^ + eg sin 0; [A — cr 



Now we determine the total effective drift velocity in an 
axisymmetric case: 



[V('')xB] 
where 



= [V^^^xBr+V^-^^xBp], , (37) 



Therefore, the effective drift velocities, V^^^ and V("*\ 
for the toroidal and poloidal magnetic fields are different. 
The additional effective drift velocity, V'^'^^ , is a result of 
an interaction of turbulent convection with inertial waves 
and Rossby waves. Indeed, a part of the tensor a[^^(k) cx 
V'nV'fiJ where -i/jj^ — 2(J7 • k)/fc is the frequency of the 
inertial waves and -i/;^ — 2AQ,xky/k^ is the frequency of 
Rossby waves [see Eqs. and ( |A37| )] . 



IV. DISCUSSION 

In this paper we studied an effect of rotation on a de- 
veloped turbulent stratified convection. This allowed us 
to determine the dependencies of the hydrodynamic he- 
licity, the alpha-tensor and the effective drift velocity of 
the mean magnetic field on the rate of rotation and an 
anisotropy of turbulence. We demonstrated that in a 
turbulent convection the alpha-effect can change its sign 
depending on the rate of rotation and an anisotropy of 
turbulence. We found different properties of the effective 
drift velocity of the mean magnetic field in a rotating 
turbulent convection. In particular, a poloidal effective 
drift velocity can be diamagnetic or paramagnetic de- 
pending on the rate of rotation. There is a difference in 
the effective drift velocities for the toroidal and poloidal 
magnetic fields which increases with the rate of rotation. 
We found also a toroidal effective drift velocity which can 
play a role of an additional differential rotation. 

Some of the results obtained in our paper using the r- 
approximation are observed in the direct numerical sim- 
ulations of the stratified turbulent convection (see p3[). 
In particular, it was found in [p3| that the alpha-effect 
can change its sign depending on the rate of rotation. It 
was also demonstrated in |^ that there is a difference in 
the effective drift velocities for the toroidal and poloidal 
magnetic fields, and that an observed toroidal effective 
drift velocity in |^ can play a role of an additional dif- 
ferential rotation. 

Now we apply the obtained results for the analysis of 
an axisymmetric afJ-dynamo. The mean magnetic field 
in an axisymmetric case is given by B = _Be^+ Vx (Ae^), 
where A is the vector potential. The equations for B and 
A in dimensionless form are given by 



dB 
'dt 
dA 



+ r±V-{Y^^''> rl^ B) ^ D {(lA) + A^B , (39) 
rji (V('4).V) (r_L^) = aB + A,A, (40) 
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where the length is measured in units of the thickness 
of the convective zone Lj,, the time is measured in units 
of L'j^/rj^, the velocity is measured in units of rj^lLc, 
the turbulent magnetic diffusion 77^ — IqUq/S, and uq 
is the characteristic turbulent velocity in the scale Zq, 
D — RaRuj is the dynamo number, Ra = Lca*/'7r 
and Ru, — Ll{5Vl)<,/ri^. Here a is measured in units 
of the maximum value a* of the a effect, ((Jfi)* is the 
characteristic differential rotation in the scale Lc, — 
[V(5rj)x V(r_L A)] • e^, = A - r'^ , r± = rsm9 and 
we used the induction equation for the mean magnetic 
field (see, e.g., |l |, §, |, §) and Eqs. d^and Q. 
When Vf-^) = V^) and V-V^-^) = 0, Eqs. (|9|) and^) 
coincide with that given in 1^ . Now we seek for a solution 
of Eqs. ( |39| ) and ( ^ ) in the form A, B oc exp(7i + ik-x), 
where k = kek , = x e^^ , the unit vector e^j is di- 
rected opposite to V(^ri) and 



±[(K/2 + ifcC/(3))2 +ifci:)]l/2 



(41) 



K = -V-V(^), f7(i'3) = V(i-3).efc and ^ ^ ^^+iuj^. 
In the limit of large dynamo number \D\ the maximum 
growth rate of the mean magnetic field 7^^ is given by 



7, = (3/4)(|i?|/4)2/3 + «,/2. 



(42) 



which is achieved at the wave number km = 
(l/2)(|D|/4)^/'^. At this wave number the frequency uj^ 
of the dynamo wave is 

= -m/^r/' - (l/2)C/(i)(|i?|/4)V3 (43) 

(see [^). The negative sign of implies that the 
dynamo waves propagate to the equator in agreement 
with the solar magnetic field observations. On the other 
hand, the divergence of the effective drift velocity V*^^) 
of the toroidal magnetic field can cause an increase of 
the growth rate of the mean magnetic field when k > 0. 
The change of the sign of the a-effect depending on the 
rate of rotation and anisotropy of turbulent convection 
(see Section III-B) can explain the observed direction of 
propagation of the solar dynamo waves. 

Note that a meridional circulation in the solar con- 
vective zone can also cause an equatorward drift of the 
solar dynamo wave (see, e.g., |43|). However, it was 

shown recently in 144] that the meridional velocity, which 
is required for the equatorward propagation of the solar 
dynamo wave with the period ^ 22 years, should be of 
the order of '--^ 10 — 12 m/s. Such large meridional veloc- 
ities are not observed on the solar surface. On the other 
hand, we found that the effective drift velocities of the 
mean magnetic field have a meridional component (along 

Be). This velocity has the maximum (vj ^^)„iax ~ 10 — 12 
m/s in the upper part of the solar convective zone. There- 
fore, this meridional effective drift velocity of the mean 
magnetic field can cause the equatorward propagation of 
the solar dynamo wave in the upper part of the solar 
convective zone. Note that the meridional circulations 



in the solar convection zone and the meridional compo- 
nent of the effective drift velocities of the mean magnetic 
field are different characteristics, because the first veloc- 
ity describes large-scale fluid motions (which may cause 
advcction of the mean magnetic field by the large-scale 
fiuid motions, i.e., by the mean fiow), and the second ve- 
locity determines the drift velocity of the mean magnetic 
field (which is originated from the mean electromotive 
force £ — (uxh)). 

We found also that in the upper part of the solar con- 
vective zone the a effect does not change its sign, i.e., it 
is positive. But in the lower part of the solar convective 
zone the a effect changes its sign, because the parame- 
ter riTg increases with the increase of the depth the so- 
lar convective zone, and the a effect becomes negative. 
Therefore, in the lower part of the solar convective zone 
the negative a effect is responsible for the equatorward 
propagation of the solar dynamo waves. On the other 
hand, the meridional effective drift velocity of the mean 
magnetic field in the lower part of the solar convective 
zone is very small and, thus, it cannot be used for the 
explanation of the equatorward propagation of the solar 
dynamo wave. 

Therefore, both effects, the meridional effective drift 
velocity of the mean magnetic field in the upper part of 
the solar convective zone and the sign reversal of the a 
effect in the lower part of the solar convective zone, can 
cause the equatorward propagation of the solar dynamo 
wave. 

Note that in the present study we did not discuss the 
magnetic buoyancy effects which play an important role 
in a creation of strongly inhomogeneous magnetic struc- 
tures (see, e.g., [i H || HH). 
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APPENDIX A: DERIVATIONS OF EQS. (|17|), 

(Ei), (Ei), @) AND (Esh. 



1. The conservation equations 

Equations (||) and yield the following conservation 
equations for the kinetic energy Wu = PpU^/2 and for 
Ws = p,Sy2 : 



dWu/dt + V • 
dWs/dt + V • F5 



Du , 



(Al) 
(A2) 



where the source terms in these equations are = 
— Pj,(u • g)S and Is ~ —lu^l/g^, the dissipative terms 
are £>„ = -Po(u ' f^) and Ds = PgSV ■ F„, the 



10 



fluxes are F„ = u (W„ + P) and F5 = uWs- Equa- 



tions (Al) and (A2) yield a conservation equation for 
We = WuCll/g^ + Ws 

dWE/dt + V -Fe = -De, (A3) 

where the dissipative term is De — Du^l/g^ + Ds 



and the flux is ¥e = Fu^l/g^ + ^3. Equation (|A^) 
does not have a source term and it imphes that with- 
out the dissipation [De = the value / We dV is con- 
served, where in the latter formula the integration over 
the volume is performed. For the convection 17^ < 



and, therefore, Ws ~ Wu\^'l\/g'^- Averaging Eq. ( |A1| ) 
over an ensemble of fluctuations we obtain a relationship 
between the flux of the entropy and the dissipation of 
the kinetic energy in a stationary turbulent co nvection: 
{uiS)gi = (u • f;^). Similarly, averaging Eq. (A3) over an 
ensemble of fluctuations we obtain (u-^) = {S"^) {g"^ / . 
Equation (Al) yields the relationship between $* and /*: 
/. - 2XgT,^:/e. 



2. Modification of turbulent convection by rotation 

Now we study a modification of turbulent convection 
by rotation. To this end we derive equations for the fol- 
lowing second moments: 

/u(k) = L{vi,Vj) , x(k) = L{w,v^) , 
Fik) ^ L{s,w) , G(k) = L{w,w) , 
$,(k) = Lis,v,) , e(k) =L(s,s) , 

using Eqs. (||)-(0), where L(a, b) — (a(k)6(— k)) and v = 
y/ po{z) u. The equations for these correlation functions 
are given by 



(A4) 



5x(k) 
dt 



1 



(*fc^„-7^J(/(k)--^G(k) 



-iV«x(k) + gi(k)F(-k) + x„ , (A5) 



dt 

dF(k) 
dt 

dG{k) 

dt 
dQjk) 

dt 



g 

-iV'«$.(k)+gx(k)e(k)-f $Ar, (A6) 

-^x(-k)-(*fc^, +^J<i>,(k) 
9 

+Fn , (A7) 
2/cV„x.(k)-2^,XH(k) + Gjv , (A8) 

(A9) 



-2-!i$7?(k) + e 

9 



N 



where 



^Rik) = [$,(k)-)-<i>,(-k)]/2, 
$/(k) = [$,(k)-$,(-k)]/2z, 



and similarly for other second moments, /„, Xw: ■ 
are the third moments which are given by 



e 



N 



/«(k) 


= L{Vn, 


v.) - 


f L{v2, Vn) , 


X„(k) 


= L{Wn 


,Vz) 


+ L{w,Vn) , 


$7v(k) 


= L{Sn, 




^L{s,Vn) , 


Fn{^) 


- L{Sn, 


w) -1 


-Hs,Wn) , 


Gjv(k) 


= L{Wn 


,w) 


-hL{w,WN) , 


ejv(k) 


= L{Sn, 


s) + 


Lis, Sn) , 



and 



Vn = -Vp^e.{Vx[Vx((u. V)u-f,)]},(A10) 
T^AT = v^e- [Vx(uxw + f^)] , (All) 

s 



S 



N 



-Vp^\ (u- V) 
1 



To 



■ div 



(A12) 



tlj^ = n-A, Vn = 2(J7-k)/fc, = 2ArJa;A:y/fc^ and 

5_L(k) = g{k±/kf. We assumed that {l/A)K^ < k"^ . Now 
we introduce the following variables: 

Xp(k) = fcV'nXH(k)-f V'AXi(k) , 

Xm(k) = fcV'n Xi (k) - V-A XflXk) , 

Fp(k) = fcV'n^i?(k)-V'AT/(k) , 

F„(k) = k4>,,Fi(k) + i,^FR{k) , 



which allow us to rewrite Eqs. ( A4)-(|A9|) as follows 
a/(k) 



dt 

gXp(k) 
dt 

^Xm(k) 

dt 

d^Rlk) 
dt 

d<^i{k) 

dt 
dFj,{k) 

dt 
dF„,{k) 

dt 
dG{k) 

dt 
dQ{k) 

dt 



-^X™(k) + 2gi(k)$^(k) 
V,X™(k)+gi(k)Tp(k)+x(f) 



(A13) 
(A14) 



{k^,,Y^f{k)-^G{k) 

-gdk)FM+x^:r^ , 



^T,„(k) + v„$/(k) 
Kg±(k)e(k) + $f) 



'/'„Xp(k) 

(A15) 



- Af 



-^F,(k)-V,<i>«(k) + <i>(^) , 



(A;VJ2<i>,(k)+4'^ 

-(fcvj2$«(k) + 4") 



— Xm(k) -I- Gat 



e 



N 



(A16) 
(A17) 

(A18) 

(A19) 

(A20) 

(A21) 



where we neglected small terms proportional to i^l/g. 
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Next, we use the T-approximation which allows us to 
expre ss the third moments /„, xjf \ • ■ • j ©Af in Eqs. 
( A13| )-( A21) in terms of the second moments [see Eqs. 
(B|)], where the superscript (0) corresponds to the back- 
ground turbulent convection (it is a turbulent convection 
without rotation, O = 0), and T(fc) is the characteristic 
relaxation time of the statistical moments. We consider 
the background turbulent convection with x^"^ (k) = 0. 

We assume that the characteristic times of variation of 
the second moments /(k), Xp(k), ■ • ■ , 6(k) are substan- 
tially larger than the correlation time T{k) for all turbu- 
lence sc ales. T his a llows us to get a stationary solution 
of Eqs. (|AT3|)-(|A21|): 



/(k) 

Xjk) 
X.(k) 

$fl(k) 

*/(k) 

i^fl(k) 
Fi{k) 
G(k) 

e(k) 



/(°)(k)-2^^[^^(k) 
+r(fc).g^(k)$fl(k)] , 
-V^AMi(k) + kip^'ip^fi20^) 

fcV'nA'lCk) , 

$W(k) 

i^f (k)-fc^,$fl(k), 

G(")(k)+2(fcVjVi(k), 

e(°)(k). 



(A22) 
(A23) 
(A24) 

(A25) 

(A26) 

(A27) 
(A28) 
(A29) 
(A30) 



where we changed t^j^ ■0^, t^j^ i/j^, Tip^ 

/xi(k) = -^^[£/(")(k) 

-T(%^(k)$^(k)(l-2V.^)] , 



Ai2(k) = Aii(k) - 



r(fc).g^(k)a>fl(k) 
1 + 



and/(°)(k)-G("nk)/fc^ = -ef^°'>(k). Here we neglected 



the terms ~ OUMq) 
term in Eq. (A23|), 



']. We will show below that the first 
X^-^)(k) = — '0^/xi(k), contributes 



to the Qf-effect, whereas the second term in Eq. (A23), 
X^^''(k) = kip^ip ^^2{^i contributes to the additional ef- 



fective drift velocity. Thus, Eqs. (A22)-(A30) describe a 
modification of turbulent convection by rotation. 



3. The correlation tensor of velocity field 

The functions /(k), G(k) and x(k) determine the cor- 
relation tensor /y (k) = {vi{k.)vj{—'k)): 



,/J)(k) + /JHk) 



(A31) 



/(k)--G(k))i^/-^(fc^) 



/J^(k) 



+ {i/2e){k,K, - fc.A,)/(k)| , (A32) 

(l/2fci){[*(k.e)i?|f)-(A/2)i?lf 
-z2fcie,,pfcp]x«(k) + [(k.e)4f) 



-(zA/2)i3(f ) - 2fciA^f^]x,(k)} , (A33) 



and Pij{k) — 
k^ + (k • e)e, 
(ki),(ki),/fci, 
(k^xe),e„ b\P 



I , and j4 



- 

1^1] - 



e^e 



(k^xe)j(kx)j 



(kj_xe)iej + 
(k^xe)j(k^)i and 



(M) _ 



= (k^xe)i(k^)j — (kj_xe)j(kx)i. For the deriva- 



B,, _ 

tion of Eqs. (A32) and ( |A33| ) the velocity is writ- 
ten as a sum of the vortical and the potential compo- 
nents, i.e., \'±_ ~ Vx(Ge) + V^(^, where v — vj_ + u^e, 
w = -A^G, Aj_if = - dvjdz, V • v = 

(A/2)(ve), V±_ = V-e(e-V). We also used the iden- 
tities (k_LXe)i(kiXe)j = k\p'>^\ki 
k'^j^^j ^ = k'j^Sijpkp (see, G.g., 
neglected the terms ~ 0[(AZo)^] 



and (k-e)B^f^ 



In Eq. (|A33| ) we 
We win use Eqs. ( [A3^ ) 
and ( [A33| ) for the calculation of the hydrodynamic hclic- 
ity and the a effect. 



4. The hydrodynamic helicity 

Now we find the dependence of the hydrodynamic he- 
licity x^"' = (u • (Vxu)) on the rate of rotation and 
anisotropy of turbulence. In k-space the hydrodynamic 
helicity is given by 



X 



(v) 



(k) = -iei„mfci/m„(k)exp(Az) 

= (l + fcVfci)xJk)exp(Az) , (A34) 



where we used Eqs. {]^^) and ( |A33D . The function 
exp(Az) in Eq. ( A34[ ) implies that we us ed the transfor- 
mation u = exp(A2:/2)v. Equation (A34) can be rewrit- 
ten as 



x('')(k)=exp(Az)[xi(k)-|-x2(k)] 



where 



XI (k) 

X2(k) 



(A35) 



(A36) 



-v,/^i(k)(i + fcVfci) , 

-4T2f2^A(t2) X e)ma)„fcm„/X2(k)(l 
+fcVfci) , (A37) 



where i2) — and we used the identity kip^ip^ = 

— 4t^J7^A((I> X e)m<^nkmn- The integration in k-space in 
Xi,2 = /xi,2(k) dk yields 



1 



Xi 



12(5, 



sin0,{(a-K3)0,{41} 
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+ (7a - 9)0,{Mit} - 3(^ - l)0,{e„„M(?)} 
+ (A/2)[Mii),(2a;) - 2I0^1{2lu)]} , (A38) 
X2 cx (cD X e),„(2;„[4P)(^) - M^l{u;)] = , (A39) 

where uj = 4r„ri, Uq = 2(7Tg<i>*(5*, sin^/ = a) • e, 
lo = uoT,, = 2X(2c^) - X(u;), e.g., = 



and 



(2) , 



/if^H = (e/TT^^'+i) / yPh,iy')dy, (A40) 



Mif)(c.) = (e/TT^^'+i) / z/Af,,(2/2)dy, (A41) 



(z) are deter- 



mined by Eqs. (Bl) and (P3) in Appendix B, and the 
exponent p = 1,2,3,4 is determined by in the ex- 
pressions for the the hydrodynamic hehcity, the a effect 
and the effective drift velocity (see below). For example, 
p = 1,2 in Eq. ( |A38D . Equation ( [A3^ ) yields the an- 
gular velocity dependence of the hydrodynamic hclicity 
X^'"^ which is given by Eq. (p7|). 



5. The electromotive force 

In order to derive equation for the a-tcnsor we intro- 
duce the electromotive force 

£, = (u X b), = exp(Az)e,™„ J xHi^) dk , (A42) 

where XifO^) — L{vi,hj) = (wi(k)/ij (— k)) is the cross- 
helicity tensor. Using equation for v = Vx(Ce)+Vx(p+ 
VzG we obtain 

-fc,e,)(fc,+zA/2)]xf (k)} , (A43) 

where xfH^) = e,x!j^(k) = L{vz,hj) and |(k) = 
L{w, h). Using Eqs. and ( pl| ) we derive equations 

for ^(k), x^'^Hk) and C(k) = i(s,h): 

^ - + (ifcV. - i'Jx'^"^ (k) -I- 1„ , (A44) 

= +k-\tk^,+^,)i{k)+t^,X^'^ 
+.9±(k)C + xL'^ (A45) 



dt 



(A46) 



where £ , v and C are the third moments: 

|„(k) - L(M^Ar,h) + L(7«,H^) , 

X^^\k) = ^(y^^,h)+^(«„HJv) , 
C„(k) = i(5jv,h)-f L(s,H^) , 



Hat = VPo Vx(uxb-£:-r;Vxb) , (A47) 



^j^^n- A, 



2(r2 • k) /fc, = 2An^ky/k'^, and 

dhi 

■I = -i(hJ-kJe./y(kj 

(A48) 

(A49) 



^(.""^'^J =-z(B.k)e.^(k) 
-A/(k)B,-, 

lL,^) =-t{B-k)x^-\k) 



Ax(k)B, 
A$.(k)B 



L s 



(B • k) *(k) 



(A50) 



k = k + (iA/2)e and x(™)(k) = £(w,v) is given by 

X'^^H^) = fcI'{x(k)[e(fc2-^AfcJ2)-k(fc, 

-iA/2)] - iG(k)(kxe)} . (A51) 

Note that x(k) = x''"'''(k) • e. Now we use the r- 
approximation and assume that the characteristic times 
of variation of the second moments £, C and x''^-' are 
substantially larger than the correlation time T{k) for all 
turbulence scal es. T h is all ows us to get a stationary so- 
lution of Eqs. (|A44|)-(|A46|): 



£(k) = rl(-)(k) + (tk^P,, - ^Jx(^Hk) , (A52) 
(c) ^ lO+ii+^rjW 

+k-\tk^p,,+^Jl(-'Hk)], (A53) 

and C(k) — rl'^''^ where we changed — > ■0^, 
Tip^ -0^;, Tip^ ip^. Now we take into account 
that a general form of the electromotive force is given by 

= ayS, + (V('^)xB), -7fe(VxB)J-K,,fe(9B)y■ 
(A54) 



-[(5x(VxB)],; = QijBj + bijkBi^j , 



(see, e.g., |4^), where the tensors and r]ij describe 
the a-effect and turbulent magnetic diffusion, respec- 
tively, V^'*' is the effective diamagnetic (or paramag- 
netic) velocity, Ktjk and S describe a nontrivial behavior 
of the mean magnetic field in an anisotropic turbulence, 
= \7jB„ (9B)„ = (l/2)(B,j ,). In this study 

we determine only the tensor cty and the velocity V^'^\ 
The calculations of the other coefficients defining electro- 
motive force is a subject of a separate paper. The tensor 



- (S) 

— '^ij 
where 



follows from Eqs. (|A43|), (|A5^)-(|A54|) 



(AS) 



''ij 
(d) 



(fc) 
(e) 



(A55) 
(A56) 
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where a[f'' = {l/2){aij + aji) and a^^'''' = {1/2) {aij -a ji) 
are the symmetric and antisymmetric parts of the tensor 
Uij, and 



aL.^(k) 



= %VASi(k,z)[2Aii(k)+/(k)-G(k)/fc2 

+2T(fc)gx(k)$fl(k)] , (A57) 
= {cikj + ejfci)s2(k, z)[V'nG(k) - fcXj(k) 
+fcT(fc)gx(k)i^/(k)] . (A58) 
-/ /(k) + G{k)/e + 2r(fc).g^(k)$H(k) 
A 1 + 



-2Ai2(k))fcV. -X.(k) 



«!f(k) 



al;^(k) 



+ [{exk),kj + (exk)jfc,]s2(k, z)[fc/(k) 
+fcT(fc)g^(k)$;^,(k) - ^.X.(k)] , (A59) 

[2fc SijmGnPmn{k) + (exk)jfcj 

-(exk)jfcj]s2(k, z)[fc/(k) 

-VnX.(k) + fcT(fc)gi(k)$fl(k)] , (A60) 

+e,fc,]s2(k,z){fc7/;J/(k) 
+r(fc).g^(k)$;,(k)]+x,(k)}, 



and si(k, z) = exp(Az)r(fc)(fc//c^)V(l + V'^), 



(A61) 
S2(k, = 



(A/2fc'^)si(k, z). Here we used that 

S = fcj2[iCx(exk) + (exx^'=^)(fc^ +iAfcj2) 
-(kxx("))(fc,+zA/2)] , (A62) 
e,/y(k) = k-^{[k^e,P,,{k)+i{A/2)hP,jie)]f{k) 

+z(exk),[x«(k)-zx,(k)]}, (A63) 



where Pij{e) — Si 



j. Note that ei/y (k) 7^ ei/ji(k) 



because rotation causes a nonzero hehcity in the turbu- 
lent convection. Here we also took into account that the 
tensor Uij must be real in r-space. 

We will show that the tensors a^J-* (k) and a^^"* (k) con- 
tribute to the a-tensor, the tensor a\'^\'k) contributes to 
the effective drift velocity V^^^ the tensor a|j^(k) con- 
tributes to the effective drift velocity Y^^") , and the tensor 
a-^'(k) contributes to the effective drift velocity V^'^^. 



The a-tensor 



Now we determine the tensor ai 
integration in k-space yields 



1 



65* \ Lp 



+3(f7 - l)0,{Af^f }] sin^bi , (A64) 



1 rim 



6S^ V L 



" ijmn IrG 



+3(a-l)03{M(f)}] 



(A65) 



where Eqs. ( |A57| ) and ( [A5^ ) for a|°^(k) and a[f(k), re- 
spectively, Pl;jl^n = Cjm{ei5nj + ejSni), and hereafter we 
use the following functions: 



2A:(2w) - X{uj) 



n 



AX{2uj) - -X{uj^) , 

TT 

-X{lo'^) - 2X{2uj) - X{lo) 



K{x} 
K{x} 



TT t: \ da 

(3-a)0,{X(3)}-(A/2)[4X(2)(2c.) 

(3-a)03{X(3)} + (A/2)[2X(2)(2u.) 

(3-(7)0,{x(2)}- A[X(i)(2cj) 
~(1 + £-1)xWH], 
(3-a)0jX(3)}_(A/2)[0,{X(2)} 

4AX(3)(2w) - (A + 25*)x(3)(a,) 



65, 



For example, 

'/-slMi^ll = -Af™„(^2) - 2M(f) (2^) - M(?) H , 

TT 

</'eU™} = (3-a)4{/(?)} + (A/2)[2/(?)(2a;) 



the f unctio ns Im ni'-^) and Mmn{^) are determined by 
Eqs. ( |A4(]|) and ( |A4lD , Mij{y) = emnlijmniy): and the 



(3), 



fimctions /ij (z) and/ym„(z) are determined by Eqs. (Bl) 
and (B2) in Appendix B. Now we use the following iden- 
tities 

P^^L^Mran = (e^tij + e-,c2>i) (Z3 -f Z2 sin^ 0/ ) 



ijmn 



+4: e^j L3 sin 



where Lfc are determined by Eqs. (C3) in Appendix C. 
Thus, the a-tensor, aij = o-J^ +a-^\ is given by Eq. (|2^). 
For a slow rotation (w <C I) the tensor Uij is given by 



a. 



5 V Lo 



— [l + (2e) 1] ]ie,^,+e,u;i) 



(A66) 



and for ^ I it is given by 



16 Up y 



ujij sin A + 
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+2cos^ 



(A67) 



a[fBj = [V(3) X (Bp - Bt)],, where B = Bt + Bp is the 
mean magnetic field with the toroidal (Bt) and poloidal 
(Bp) components, the tensor a[f {h) is determined by 



Eq. ( A5S| ). Integration in k-space, a\f — J a,-^^ (k) dk, 
yields 



7. The effective drift velocity 



Now we determine the effective drift velocity V, 



id) 



Vr 



(1) 



and 



Vm \ where 



1/(1) 



-(l/2)e„.,a(f = 
-(l/2)e™,a(f , 



pid) ~{d) 
mij ij 



(A68) 
(A69) 



4f = -^(^^^{2cj2(u;xe)„c2)„[0J/.yw} 

+ 2(3 - (T)(f)^{Iij„in} + 6{a - l)(j)^{Jijmn}] 



where P^/^in = i^ipnSjm + £jpnStm)ep. In order to deter- 
mine the effective drift velocity V^'^' we use the following 
identities: 



Ad) 



~(d) 

n 



-Op.. A \7,id:) 
1 / ZoWo 



+3((7-l)0,{M^/a] , 

1 
6^ 

+3(a-l)0,{Mi3)}] 



(A70) 



(A71) 



(A72) 



P. 



id) 



ej5in)'jJm- For the integration in k-space we used 
Eqs. ( |A60| ) and ( |A61| ) for a'^\\<i) and af^{\<i), respec- 
tively. Using the following identities: 



hid) J 
p(d) iTr 



= -CiLi + LJi sin A2 , 



e,[3Ci-Ai-f (3C'3-A2)sin2 0,] 
+tDi sin (j>i [36*3 -I- C2 siv? , 
Im.n. = -2Zi(d>xe)fc , 

-2(L3 + L2 sin^ (/)i)(tDxe)fc 



ties f/^'' and V^^' which are given by Eqs. (^5|) and (26). 
Here Mij{y) = e„^n£^^^mr^(?/), ^(2^) and hjmniz) are de- 
termined by Eqs. ( |BlD a nd ( |B2[ ) in Appendix B, ifc are 
determined by Eqs. ( |C3| ) in Appendix C. 

The electromotive force has a term Bj which for an 
axisymmetric case contributes only to an additional effec- 
tive drift velocity, V(^\ of the mean magnetic field, i.e., 



in Eqs. (A68)-(A72), we obtain the effective drift veloci- 



iqj+c.jqi)Bj = [(qxc)x(Bp - Bt)], , 

Pijlin^mn = -(^3 + C2 Sin^ (l3l){CiLUj + CjLOi) 

-2(73 sin (c^e-, +Cjei) , 



(L5 + Lq sin^ (j)i){ctLjj + CjLUi) 
+Di sin 0; (cjCj -I- CjCj) , 

where = (a>xe)i, q = a> or q = e, Afij = Jijmm, 



rip) 



{u:) = (6/^c^P+i) / j/%>„„(y^)dy , (A73) 



= {^1^^"^') / y^4™„(y^)dy,(A74) 



and we used Eqs. (B1)-(B9) in Appendix B and Eqs. 
(|C2|)-(C5 ) in Appendix C. Thus, the effective drift veloc- 
ity V'^' is given by Eq. (^. 



APPENDIX B: THE IDENTITIES USED FOR 
THE INTEGRATION IN k~SPACE 

To integrate over the angles in k-space we used the 
following identities: 



hj (a) 



, ^^'"'''2. de = Ai(5,, + A2 oj,, , (Bl) 
1 + a cos^ 9 



kijnm sin 

I + a cos^ 9 

+C3{dijUJmn + Sim^jn + SinUJjm + Sjmi^in + ^jn^im + Smn^ij) , (B2) 



lijmnia) — I tt-t: d9 dip ~ Ci{SijSmn + SimSjn + Si.n6j,n) + C2l^ij 

' 1 + a cos^ '3 
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Jijmn ('^) 



J\i ' ' Sin — 1 — 

epq I -['^™^^^Qg2 '^(^dip = [Di + -Dsiu; ■ e)\SijSmn + StrnSjn + StnSjm) 
+ [D2 + Dt{u} ■ ef]uJijmn + D^{u} ■ e){uJtj,nen + (^ijnem + i^jmnei + LOimnej) 



^iimn sin ^ 

-— ^— dfc' dip = 

(l + acos2 6')2 ^ 



do J + cos^ J 



5 

%„„sin6' - ad- 

— — ; „ «P a(y5 = -tlijmn[a) + 7; J^ijmn(a) 



My (a) = (Ci + C3 sin^ + (C3 + C2 sin^ 0;) + 2Ci e. 



+2C3 sin {ciUOj + gj Wi) 



(B3) 

(B4) 
(B5) 

(B6) 



where My (a) = eTO„/y„m(a), 

e.jMy (a) = SCi + 66*3 sin^ 0, + C2 sin"* , (B7) 

Mpp{a) EE 6^4- (a) = Ai + A2 sin^ cj)i , (B8) 

and LOij = uJiUij, uj^jm = LUiLUjLUj, Ai = 5Ci + C3, 
A2 = C2+7C3, and 

Alia) = F(l;-l;0;0),i2(a) = i^(-l;3;0;0), 

Ci(a) = (1/4)F(1;-2;1;0), 

(72(a) = {l/4)F(3;-30;35;0), 

C3(a) = (l/4)F(-l;6;-5;0), 

Ci{a) = i^(0;0;l;-l) , (75(a) = F(0;0;-1;3) , 

Di{a) EE -(l/8)((7i + 5(73 - 5(74) 

= (l/8)F(l;-7;ll;-5), 
D2(a) = -(l/8)(51(7i + 111(73 - 119(74) 

= (1/8)F(15;-141;245;-119) , 
^3 (a) = (3/8)(3(7i + 7^3-7(74) 

= (3/8)F(-l;9;-15;7) , 
^4(a) = (l/2)((7i + C3-(74) 

= (1/2)F(0;1;-2;1), 
55 (a) = 3(7i + 9C3 - 7(74 

= (l/2)F(-3;24;-35;14), 
De{a) EE (l/8)(5(72 + 3(73 + 20(74 - 5(75) 

= (l/8)F(3;-33;65;-35) , 
57(a) = -(l/8)(48(7i + 27^2 + 165(73 + 28(74 

-35(75) = (1/8)F(9;-21;-105;133) .(B9) 

Here 

F{a; P; 7; fl) = TT[aJoia) + /3J2(a) + 714(0) + ^Jeia)] , 



J2fc(a) = 2 / a;2V(l + aa;2)t^a; = a~i[2/(2/c-l) 
Jo 

-J2(fe-i)(a)] , (BIO) 



and Jo(a) = 2arctan(Y^) / ^/a- In the case of a ^ 1 
these functions are given by 



J2k{a) 



2k + 1 



1-a 



2k + 1 
2fc + 3 



,2fc + 1 
2fc + 5 



and for a ^ 1 they are given by J2fe(a) ^ 2/a {2k — 1) for 
ah integer k except for fc = and Jo(a) ^ T:/^/a — 2/a. 
Now we introduce the following functions: 



where 



EE + p4^\u) + ^j'i'\u^) + 



41\lo) = (e/c^^'+i) / yKhuhi") dy . (Bll) 



The integration in Eq. (Bll) yields: 
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(2fc-l)(p-l) 



J^f;_^j,(c.) (B12) 



for p 7^ 1 and all integer fc except for fc = 0. When p = 1 
and fc ^ we get: 



6 



2fc- 1 



ln(l +^2) 



E 



m, ,2m 



2 / arctan(a;) 



(-l)™a; 



m=0 

When fc = we obtain 

6 

n 



(2m + 1) 



(B13) 



arctan(a;) / (-l)"+i 

^ -|^^n+m— 1 



n-1 



(2m- l)w"-™+i 



(B14) 



16 



J, 



(2n+l) 







6 



2n+l 

ri-l 



,arctan(a;) , (-l)"+i ^ 
■ ln(l + to ) 



-n 



Cj2("+1) 
,„ (1+Cj2)»-™_1 



m— 1 



(n — m) m! (n — m)\ 



(B15) 



Equation ( [B14[) is for all integer n > 1, and Jq^\uj) = 
(12/ix)) [arctan(y)/j/] dy. For n = and n — 1 the third 
term with the su m in Eg. ( |B15| ) should be dropped. In 
order to use Eq. ( |B12 ) p — 2 we need to know the 



function J2^f} (lo) which is given by 



2k-2 

fe-l 



arctan((jj) f ^ ^ (— 1)*^+^ 



,2k 



(-1) 



^ (2m + l)w2(fc-m) 

In the case of a; ^ 1 these functions are given by 



(B16) 



12 


1 - 


2 

- LO 


(2fc4 


^1)(P+1) 


+L0^ 


(2k + l\ / 






{2k + 5)\ 




15) 



f2k + l\fp+l 



\2k + 3j \p + 3 



In the case of a; ^ 1 these functions are given by 



(2/fc- l)(p- l)cj2 



for p ^ I and k ^ 0; 



for p and p ^ 1; 



12 



puj {p — l)uj^ 



for fc 7^ 0; and 



12 In ( 



4'\iu) 



2*^-2 
67r 



(2fc- l)w2 
1 /4fc-l 

TTCJ 



2 In u) 
1 



2k- 1 



67r In UJ 



Now we introduce the following functions 



(p) 

ijmn 



(co) 



G^"^ (uj) 



Jo 

(3/7r)4„„,(a,2) - {p-l)llf^Jco)/2 , (B17) 

(e/TTL^f+l) r yPG,,mn{y^)dy 



P_Z1\ T^P) 

2 / ijmn 



3cij2 ^ dIijmn{o) 
TT 



da 



(B18) 



which will be used for the calculation of the effective 
drift velocity of the mean magnetic field. The functions 
i(p) 

of LHS of Eqs 



{uj) can be obtained from Eqs. (B9) after the change 
^) Ak{a) A^j^\uo) and of RHS of 
Eqs. ® F{a;P;f,fi) ^ /3; 7; /i) and similarly 

for the functions Cl!'\uj) and d'^^\ljj), e.g., 



A'i\oj) = f(^')(1;-1;0;0), 
c'f\Lo) = (1/4)f(?')(1;-2;1;0) 



(B19) 



and similarly for the other functions C^!'^ (lu) and i?^^'' (w). 
For the calculation of the functions we need to use 

the following identities: 



, / dJ2k{a) 



= - J2k{a) 



dJ2(k~i){a) 



\ 9a J^^ 
where 

dJoja) 
da 2a V a + 1 



da 



Jo(a) 



APPENDIX C: THE FUNCTIONS ^0{iu) AND 

The functions 'l'fc(w) are given by 

*iH = lOa0,{4')} + (a + 3)0,{4')} 
^i\/2)[5A[^\2Lu) + A^2'-\2uj)] 
-9(a-l)0jCf)}, 
(7a-9)0j4^)}-18(a-l)</),{cf } 
+ (A/2)4^)(2^) , 



^4(0;) 
*9(w) 

*io(w) 



-3(a-l)0,{cf }, 
0,{Ai} + 3(a-l)0,{cf)}, 
3(a-l)0,{cf }, 
0,{A2} + 3(a-l)0,{cf }, 
3(a-l)0,{cf }, 
6(a-l)[<^,{cf)} + 203{4^)}] 
0jLi} + 3(a-l)03{Lf}, 
2*5H + 3(a-l)03{Lf }. 



(CI) 



The functions Ek{w) are given by 

3(a-l)0,{if^}~0,{i4}, 



£'i(w 

£;2('^ 



= 0,{A2} + (3a-l)0,{if^}, 

= 3(a-l)0,{cf }, 

= 0,{A2} + 9(a-l)0,{cf }, 

= 0jLi} + 3(a-l)0,{Lf }, 
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Eio{uj) 
where 



3(a-l)0,{Lf }, 

(1/2)[<^,{A2} + 3(a {Cf}] 

W[0jL3} + 2(3-(7)0,{i3} 

+6(a - 1)(/)JL5}] , 
-Ejiuj) - Egiuj) , 

-(3/2)(a - 1)0, {Cf} - 6a>2(a - Ij^jLg} 



= £;iH + 3(a-l)[0,{cf } 



(C2) 



Li(a) EE Ai+A2 = 2^(0;l;0;0), 

L2(a) = (72 + 3(73 = F(0; -3; 5; 0) , 

Z3(a) = (7i+C3 = Z(0;l;-l;0), 

Z4(a) = 2A1+A2 =F(1;1;0;0) , 

Lsia) = ^i+^3 = (l/4)F(-l;10;-17;8), 

Z6(a) = (1/3)53 + 56 - (l/4)i^(l; -12; 25; 28), 

L^ia) = 3Ci- Ai = (l/4)Z(-l;-2;3;0) , 

Ls{a) = 6C3-A2 = (l/2)^(-l;12;-15;0), (C3) 

and Li^^u) = (6/^a;P+i) 2/^1^(2/2) dy. 



APPENDIX D: THE MODEL OF THE 
BACKGROUND TURBULENT CONVECTION 

A simple approximate model for the three-dimensional 
isotropic Navier-Stokes tm-bulence is described by a 
two-point correlation function of the velocity field 
fij{t,x,y) — {ui{t,x)uj{t,y)) with the Kolmogorov spec- 
trum W{k) oc k~'' and q = 5/3. The turbulent con- 
vection is determined not only by the turbulent ve- 
locity field u(t,x) but the fluctuations of the entropy 
s{t,x.). This implies that for the description of the tur- 
bulent convection one needs additional correlation func- 
tions, e.g., the turbulent flux of entropy $i(t,x, y) — 
{s{t,x)ui{t,y)) and the second moment of the entropy 
fluctuations 0(t,x, y) — {s{t,x.)s{t,y)). Note also that 
the turbulent convection is anisotropic. 

Now we derive Eqs. (||) and (|l^) for the correlation 
functions fij and To this end, the velocity uj^ is 
written as a sum of the vortical and the potential compo- 
nents, i.e., uj_ = Vx((7e)-|- Vj.0, where w = (Vxu)^ = 
-A±C, Aj_4> = Au^-du^/dz, Vj_ = V-e(e- V). Thus, 
in k-space the velocity u is given by 



where we neglected terms ^ 0(A). Multiplying Eq. (Dl) 
for Ui(ki) by Uj (k2) and averaging over turbulent veloc- 
ity field we obtain 



/^)(k) = fc;4[fc4/(")(k)e„„P,„(k)F,„(k) 
+(exk),(exk),G('')(k)] , 



where we assumed the turbulent velocity field in the 
background turbulent convection is non-helical. Now we 
use an identity 

{k/k±^ emnPimO^)PjnO^) = 6ij + ^ij ^ kij 

= P.,(k)-/^^(k^), (D3) 

which can be derived from 



kz{kzG-ij ^ikj -\- Gjk^ ) — k-ijk 



2 k^k"^ 



Here we also used the identity (k^xe)i(kxxe)j 



klP}^\kj_). Substituting Eq. (|d|) into Eq. (|d|) wc ob- 



^± ij 

tain 



fl°\-k) =(fc/fcx)^{/(°nk)P,:,(k) 

+ [G(")(k)/fc2-/(")(k)]/^^(k^)}, (D4) 

Thus two independent functions determine the correla- 
tion function of the turbulent velocity field. In isotropic 
three-dimensional turbulent flow G^°\'k)/k^ = /("^(k) 
and the correlation function reads 



fg\k)^ UW{k)P,,{k)/8TTk^ 



(D5) 



In isotropic two-dimensional turbulent flow 
G(OHk)/fc^ > /*/*°)(k) and the correlation func- 
tion is given by 



f^'\k)^G^°Hk)PMk^)/8nk'kj 



(D6) 



A simplest generalization of these correlation functions is 
an assumption that G(°)(k)/[fc2/^°Hk)] - 1 = e = const 
and thus the correlation function f^p(k) is given by 
Eq. (p|). This correlation function can be considered as a 



combination of Eqs. (D5) and (D6) for three-dimensional 
and two-dimensional turbulence. When e depends on the 
wave vector k, the correlation function /-^^■'(k) is deter- 
mined by two spectrum functions. 

Now we derive Eq. ( \U\) for the turbulent flux of en- 
tropy. Multiplying Eq. (pl|) written for Ui(k2) by s(ki) 
and averaging over turbulent velocity field we obtain 
Eq. iM). Multiplying Eq. iM) by j(kj_xe), we get 



*f(k) 



(D7) 



Now we assume that *i°^(k) oc *l/(°)(k)//,. The inte- 
gration in k-space in Eq. (|D7[ ) yields the numerical factor 
in Eq. (|l^) . Note that for simplicity we assumed that the 
correlation functions i^(°)(k) and /(°)(k) have the same 
spectrum. If these functions have different spectra, it 
results only in a different magnitude of a numerical coef- 
ficient in Eq. (|l|). 

Now let us discuss the physical meaning of the param- 
eter (T. To this end we will derive the equation for the 
two-point correlation function $i°-'(r) = (s(x)u(x -f- r)) 
of the turbulent flux of entropy for the background turbu- 
lent convection (which corresponds to Eq. ( pi] ) written in 



(D2) 
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k-space). To this end we rewrite Eq. (11 ) in the following 
form: 

<I>i")(k) = $*[A:2+r(e-k)2]$^(fc) , (D8) 
= -(3- CT)l¥(fc)/87rfc'' , (D9) 

where $* = **• e, T = 3(cr - l)/(3 - a). The Fourier 



transformation of Eq. (D8) yields 

$(")(r) - $:[A + r(e. V)2]$,„(r) , (DIO) 

where (r) is the Fourier transformation of the function 
Now we use the identity 



(Dll) 



where '0(r) = r J_$L(r) and ^A'C*") — dip /dr. Equa- 
tions (DIO) and (Dll) yield the two-point correlation 
functional"-' (r) : 

ln\, . N , 1 -I- r COS^ fi' \ , 

$W(r) = ^:[ij{r)+riP'{r) j ,(D12) 

where 9 is the angle between e and r. The function 
ipir) has the following properties: ip{r = 0) = 1 and 
{rip')r=o = 0, e.g., the function ilj{r) — 1 — {r/loY^^ sat- 
isfies the above properties, where 1 < g < 3. Thus, the 



two-point correlation function $1°^ (r) of the flux of en- 
tropy for the background turbulent convection is given 
by 



<i>i")(r)=<f>: 



1 - 



[q- 1)(1 + rcos2 



3-t-r 



where 1 < g < 3. The simple analysis shows that 
—il{q — 1) < (7 < 3, where we took into account that 
S^i"' (r) /dr < for all angles 0. The parameter a can be 
presented in the form = [1 -|-|(g-t- l)/(g- 1)]/(1 -t-|/3), 
where ^ = {l^/lzY^^ — 1, li_ and 1^ are the horizontal 
{9 = 7r/2) and vertical (6* = 0) scales in which the cor- 
relation function <^'f^ (r) tends to zero. The parameter ^ 
describes the degree of thermal anisotropy. In particular, 
when l±_ = Iz the parameter ^ = and a = \. For <C Iz 
the parameter ^ = — 1 and a = — 3/(g— 1). The maximum 
value Cmax of the parameter ^ is given by fmax = 9^1 
for (7 = 3. Thus, for a < 1 the thermal structures have 
the form of column or thermal jets < Iz), and a > 1 
there exist the "pancake" thermal structures {l± > Zz) in 
the background turbulent convection. 
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